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Abstract 

A difficulty in the classical hydro-dynamic analysis of moving contact-line problems, 
associated with the no-slip wall boundary condition resulting in an unbalanced 
divergence of the viscous stresses, is reexamined with a smoothed, finite-width in- 
terface model. The analysis in the sharp-interface limit shows that the singularity 
of the viscous stress can be balanced by another singularity of the unbalanced sur- 
face stress. The dynamic contact angle is determined by surface tension, viscosity, 
contact-line velocity and a single non-dimensional parameter reflecting the length- 
scale ratio between interface width and the thickness of the first molecule layer at 
the wall surface. The widely used Navier boundary condition and Cox's hypothesis 
are also derived following the same procedure by permitting finite-wall slip. 
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1 Introduction 



Immiscible two-phase flows with moving contact lines occur in a variety of 
applications, such as coating and biological processes. The moving contact 
line problem, however, has for many years remained a partially open issue. 
One of the problems is the validity of the no-slip wall boundary condition, 
which arises with classical hydrodynamics, where for a no-slip wall an unbal- 
anced divergence of the viscous stress occurs, which leads to a violation of the 
contact-angle condition at a moving contact line [8 J [5 J . There have been many 
attempts to resolve the problem by modifying the boundary condition, in- 
cluding the slip model [6j[18], the interface relaxation model [13] ; the diffusive 
interface model [10] [3], and the combined molecular-dynamics and diffusive- 
interface model [H]. However, studies by molecular dynamics show that even 
though considerable contact-line velocities can be obtained [TT] [IB], the maxi- 
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mum shear rate is still many orders less than taht which can violate the no-slip 
wall boundary condition considerably [17] . 



2 Smoothed, finite-width interface model for moving contact line 



In this Letter, we reexamine the hydrodynamics of a fluid/fluid/solid system 
with a steady moving contact line. Instead of considering sharp interfaces 
directly, our analysis starts from smoothed, finite-width interfaces. Given the 
continuous interfacial free energy density with the form [9] / = |er|VC| 2 + 
\P(C), where C is a color function, a is a coefficient and \P(C) is the bulk 
energy density, at the state which minimizes T = J fdV the interface reaches 
its equilibrium profile. In this case the surface stress Uij in a two dimensional 
Cartesian coordinate system is given by 

/ dC dC dC dC\ 

n * = T^^-^J' *>3> k = 1 > 2 c 1 ) 

where 5y is the Kronecker delta. One important property of the surface stress 
is that one of the principle axes x\ is aligned with the gradient of the color 
function, the other principle axis x' 2 is aligned with the interface tangential 
direction, along which the only non-zero component of the surface stress is 
the positive normal stress (tension), Uyy = cr|VC| 2 . The relation between the 
surface tension 7 and a for an infinite plane interface is given by 
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We consider a steady moving contact line with dynamic contact angle a, and 
velocity U s , as shown in Fig. [TJ Around the contact line there are three phases: 
fluid 1, fluid 2 and the static wall. We define the color function as 



ctl= 1 in phase/ kJ = 1 ^ (3) 
I else 

Note that the color function is discontinuous across the interfaces. In order 
to obtain a finite, continuous surface stress, we introduce a two-dimensional 
smoothing-kernel function [13] [12] 



^(x,0 = 7^e- x2 /« 2 (4) 

in which £ is smoothing length and £ << L, L is the characteristic length 
scale of the system. W(x, £) is radially symmetric and has the properties 
/ W{x. — x', £)g?x = 1 and lim^ W(x — x', £) = <5(x — x'). After convolution 
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with the kernel function, the smoothed gradient of the color function pointing 
towards phase / at a point x' in phase k is 



VC H (x') = J C H (x)V^(x - x', Odx, I + k. (5) 

Assuming that the smoothed profile defined by VC fc, (x') is the interface profile 
corresponding to an equilibrium form of the bulk energy density, the total 
surface stress at a point in phase k is can be calculated from Eq. (pQ), by 
Iljj = J2i^k n^VC^). It is easy to verify by Eq. (j2J) that, for a infinite plane 
interface between phase k and phase /, the surface tension is 

Figure [U indicates the regions of non- vanishing VC kl for different phase pair- 
ings. Note that there are overlap regions near the contact line. 

In order to study the contact-line dynamics, as shown in Fig. [T], we define a 
small square control volume with side length 2e, e <C £, with one side on the 
wall surface so that the interface between the fluid 1 and fluid 2 goes through 
the control volume. Assuming incompressibility and straight interfaces with a 
not far from |, we obtain 

2eU l n + J Y^ 1 dx l +2eT l 11 + J r^dxt = 2eU r u + J Yi{ 1 dx 1 +2eT r ll + J r^dxi, (7) 

by considering the force balance on the control volume in tangential (wall par- 
allel) direction, where the superscripts I, w, r and / represent the left, wall, 
right and upper faces of the control volume. Iln and II21 are the tangential 
components of surface stress, m and t%\ are the tangential components of 
viscous stress. As e is small and e <C £ , the gradients of the color functions 
at a point on the face of control volumes can be approximated with the repre- 
sentative values on the contact line. If the contact line is defined to be at the 
origin of a two-dimensional polar coordinate system, the gradient is given by 

VC w (x' -> 0) = / W(x, £)dx, l + k (8) 
Jn,(6,e') 

where Q(9, 9') represents the sector between polar angles 9 and 9' in two- 
dimensional polar coordinates, and 9, 9' = 0, a, it depending on the choice of 
phase pairs. It can be readily obtained that Jq^a VVK(x, £)c£x = (sin# — 
sin0', cos6>' — cosf?)^^. With Eq. (pQ), the tangential surface-stress compo- 
nents in Eq. are U\j = - cos a) 2 , II] 2 = ^sina(cosa - 1), 
n 2 l = i^(l + cosa) 2 , n 2 } = 7r ^sin«(l + cosa), U\f = ^f, n% = 0, 
n?]" = ^ and Illi = 0, where ai w , o\ w and o\i are the coefficients between 
fluid 1 and wall, fluid 2 and wall, and fluid 1 and fluid 2. Hence, using Eq. (jBJ), 
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Eq. ((7j) becomes 

1 

= l2w + ~ln c °s « 
(9) 

where 7i TO , jx w and 712 are the surface tensions between fluid 1 and wall, fluid 
2 and wall, and fluid 1 and fluid 2, respectively. Note that, the shear rates 
on the faces of the control volume are approximated with the values on face- 
centers, and that the viscosities of fluid 1 and fluid 2 are assumed to have the 
same value /1. When the fluids are in static equilibrium the second term on 
the left-hand-side disappears, Eq. becomes 

1 , 

llw = l2w + -jl\2 cos a (10) 

where a' is the static contact angle. Eq. (TTUT) implies that the static contact 
angle is different from that obtained by Young's relation, except a! — ~. 
This is not unexpected because the current relation gives the force balance 
within the interface. Note that for a 7^ a', an unbalanced surface stress along 
the tangential direction arises in Eq. ([7]), and is balanced by the differences 
between the shear stresses. To study the details of the balance between surface 
forces and viscous forces, it is convenient to define a layer, as shown in Fig. [U 
which has a small thickness of e and a velocity U E in the center of the control 
volume. As U £ — * U s for e — > 0, one can study the force balance exactly 
at the contact line. For any location other than the contact line there is no 
unbalanced surface stress as in Eq. fl7|) and hence viscous forces are continuous. 
Therefore, it is straightforward to assume that the fluid velocity on the left and 
right faces of the control volume match continuously with the wall velocity, 
and the shear rate on the upper face of the control volume match continuously 
with that of the bulk flow. Here, three types of wall boundary conditions with 
different wall-slips, i.e. no-slip, finite-slip and free-slip, are to be considered. 



Titu + \ 77^ 



dxo 



dxo 



dxi 



3 Discussion 



If a no-slip wall boundary condition is applied, the viscous stress is calculated 
from viscosity and shear rate. As e is small, a linear approximation of the 
shear rates is sufficient, then Eq. (jUJ) can be rewritten as 

TfiU e — v / 2~7r£/x ( ^— ^ J = cTi 2 (cosa — cos a'). (11) 

where T = ^> 1 is a non-dimensional parameter. Note that the normal 

viscous stresses here cancel out because of opposite directions and same mag- 
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nitudes. For the distinguished limit of a sharp interface r 3> 1 as £, e — > 0, 
one has 

rCa = cos a — cos a' (12) 
where the capillary number is defined by Ca = /iC/ s /7i2, since ff^j) ~ j; is 

finite and U £ — > U s . Now that the same form as Eq. (fT2|) with T = -y/^f f can 
be derived if the viscosity of fluid 2 is neglected. Here, the problem discussed 
in Refs. [8] and [5] can be solved: for infinite viscous stress in Eq. in the 
limit e — > 0, Eq. (j!2p implies that there still is a contact-angle condition at 
the contact line for a non- vanishing contact-line velocity U s . The reason is 
that, for £ — > the unbalanced surface stress in Eq. (J7]) becomes infinite as 
well, and Eq. ffl2l gives the condition for which the two infinite stresses are 
in equilibrium. Note that the problem of a singular viscous force remains only 
if e — > and £ does not vanish. This limit results in U s — > 0, i.e. the contact 
line does not move. 

A straightforward interpretation of Eq. ([121) with respect to the microscopic 
length scales indicates that £ corresponds to the physical width of the inter- 
face and e to the thickness of the first molecules layer at the wall surface, 
and r is just the ratio between the two length scales. Since the thickness of 
the first molecule layer is quite close to the molecule size and the interface 
width is at least several times the molecule size, our result T ^> 1 is physi- 
cally meaningful. Eq. ffl2|) is derived from a classical hydrodynamic analysis 
in which the only considered dissipation mechanism is the viscous force. It is 
quite surprising that Eq. ( fl2|) has the same form as a linearized formulation 
of the molecular-kinetic model p] which was proposed to discard dissipation 
due to viscous flow: T'Ca = (cos a — cos a'), where V = rj/fi, r\ is the coef- 
ficient of wetting-line friction. Note that r] has the units of the viscosity /i, 
and is always much larger than /i |2J, which is in agreement with our result 
r ^> 1. It is also interesting that Eq. (1121) has the same form as the small- 
velocity-approximation relation of Shihkmuraev's interface relaxation model 
[TB] . However, Shihkmuraev obtained a zero contact-line velocity for negligible 
interface relaxation, as opposite to finite contact-line velocity obtained in our 
current analysis. 

If a finite-slip wall boundary condition is applied, the viscous stress on the 
wall for e — > is given as /3U S , where (3 is the slip-coefficient. In this case Eq. 
Q becomes 



(3U S — til— — = v2ir — (cosa — cos a'). (13) 
\dx 2 J £ 

Eq. ffT31 is in agreement with the slip-wall boundary condition of Qian, Wang 
& Sheng [TJ], which states that the wall slip is proportional to the sum of 
the viscous stress and the uncompensated Young stress. An important result 
different from that of the sharp interface model is that for a given contact-line 
velocity the dynamic contact angle is strongly affected by the shear rate of 
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the bulk flow. Note that Eq. (jl~3|) is valid only for finite-thickness interfaces 
and implies a surface-force singularity if the interface thickness tends to zero. 
To eliminate the singularity the dynamic and static contact angles should be 
equal a = a', which explains the underlying reason of Cox's hypothesis [I] for 
a macroscopic analysis stating that wall-slip is permitted and the contact angle 
is independent of contact-line velocity. A result of the surface-force balance 
with a = a' from Eq. (fl~3l) is the widely used Navier boundary condition 

U s = A (|^J , where A = fi//3 is the slip length. 

Our analysis does not allow for a free-slip boundary condition [7] because free 
slip would result in the first term on the left-hand-side of Eq. ffl~3l) to vanish, 
which may lead to an un-physical an decrease of the contact angle for an 
advancing contact line. For the sharp-interface limit with free slip there is no 
force balance, no matter whether the dynamic and static contact angles are 
the same or not. 



4 Conclusion 

To summarize, we have studied the force balance at a moving contact line with 
different boundary conditions. It is found that, for the sharp-interface limit, 
both the finite-slip and the no-slip wall boundary conditions are possible. With 
the finite-slip wall assumption, the analysis explains that the previously used 
Cox's hypothesis and the Navier boundary condition are essential for a force 
balance. However, the analysis also suggests that the no-slip wall boundary 
conditions is still valid along with a contact-angle condition, which agrees with 
several previous studies on dynamic contact angles. More importantly, since 
there is no conflict with the results of molecular dynamics simulations, the no- 
slip wall boundary condition can serve for obtaining more reliable numerical 
predictions of moving contact line problems. 
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Fig. 1. smoothed, finite width interface model for moving contact line problem 



8 



